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Abstract. The development of patch test consistent quasicontinuum energies for multi-dimensional 
crystalline solids modeled by many-body potentials remains a challenge. The original quasicontin- 
uum energy (QCE) [28] has been implemented for many-body potentials in two and three space 
dimensions, but it is not patch test consistent. We propose that by blending the atomistic and 
corresponding Cauchy-Born continuum models of QCE in an interfacial region with thickness of 
a small number k of blended atoms, a general quasicontinuum energy (BQCE) can be developed 
with the potential to significantly improve the accuracy of QCE near lattice instabilities such as 
dislocation formation and motion. 

In this paper, we give an error analysis of the blended quasicontinuum energy (BQCE) for a 
periodic one-dimensional chain of atoms with next-nearest neighbor interactions. Our analysis 
includes the optimization of the blending function for an improved convergence rate. We show that 
the £ 2 strain error for the non-blended QCE energy (QCE), which has low order 0(e 1 ^ 2 ) where e 
is the atomistic length scale [131129] . can be reduced by a factor of k 3 ^ 2 for an optimized blending 
function where k is the number of atoms in the blending region. The QCE energy has been further 
shown to suffer from a O(l) error in the critical strain at which the lattice loses stability [16]. We 
prove that the error in the critical strain of BQCE can be reduced by a factor of k 2 for an optimized 
blending function, thus demonstrating that the BQCE energy for an optimized blending function 
has the potential to give an accurate approximation of the deformation near lattice instabilities 
such as crack growth. 



1. Introduction 

Crystalline materials often have highly singular strain fields at crack tips, dislocations, and grain 
boundaries that require the accuracy of atomistic modeling only in small regions surrounding these 
defects. However, these localized defects interact through long-ranged elastic fields with a much 
larger region where the strain gradients are sufficiently small to allow accurate and efficient approx- 
imation by coarse-grained continuum finite element models. This has motivated the development 
of numerical methods that couple atomistic regions with continuum regions to compute problems 
with length scales that are sufficiently large for accurate and reliable scientific and engineering 
application 0[ini[l9H2ll[23l[Ml[M[M[M[Ml[38]. 

The quasicontinuum energy (QCE) [30] couples an atomistic model to a finite element continuum 
model based on the Cauchy-Born strain energy density. The Cauchy-Born strain energy density 
reproduces the atomistic energy density for a uniformly deformed lattice. The QCE atomistic-to- 
continuum coupling also reproduces the atomistic energy density for a uniformly deformed lattice 
simply by modifying the volume of the triangles or tetrahedra within the cut-off radius of the 



Date: May 24, 2011. 

2000 Mathematics Subject Classification. 65Z05,70C20. 

Key words and phrases, quasicontinuum, error analysis, atomistic to continuum. 

This work was supported in part by DMS-0757355, DMS-0811039, the PIRE Grant OISE-0967140, the Institute 
for Mathematics and Its Applications, and the University of Minnesota Supercomputing Institute. This work was 
also supported by the Department of Energy under Award Number de-sc0002085 



1 



ANALYSIS OF ENERGY-BASED BLENDED QUASICONTINUUM APPROXIMATIONS 



2 



boundary of the atomistic region by using the Voronoi cell [28J. The QCE energy remains the 
most popular quasicontinuum energy since the modification of the volume of the triangles near 
the interface to reproduce the atomistic energy density can be explicitly achieved and has been 
implemented in the code [27] for multi-dimensional problems with many-body potentials. 

A quasicontinuum energy is said to be patch test consistent if it reproduces the zero net forces 
at each atom for a uniformly deformed lattice. The fully atomistic and Cauchy-Born energies are 
patch test consistent by symmetry, but the symmetry is broken in the QCE atomistic-to-continuum 
interface. The nonzero forces in the QCE atomistic-to-continuum interface of a uniformly deformed 
lattice are called ghost forces [HJ[T5l[28l[35] ' 

The ghost force correction method (GFC) was developed to improve the accuracy of QCE [35] . 
and it has been shown to converge to the patch test consistent force-based quasicontinuum method 
(QCF) OH]- GFC has been implemented in the code [27] and further developed to utilize the 
efficiency of atomistic-to-continuum modeling and continuum mesh a posteriori adaptivity [Tj[2"ll28t 
[35]. 

However, it has been shown that the force-based quasicontinuum method (QCF) gives a non- 
conservative force field [9~l llll[2"8] . so the stability of a deformation that satisfies the QCF equations 
cannot be simply determined by checking whether it is a local minimum of a quasicontinuum 
energy [TJ]. The linearization of the QCF method has also been shown to be indefinite and to have 
unusual stability properties [15] which presents a further challenge to the notion of stability and to 
the development of efficient and reliable iterative methods [171 125]. 

The potential for more reliable and efficient iterative solution methods and more direct extensions 
to finite temperature dynamical methods [18] make energy-based quasicontinuum methods more 
desirable than force-based quasicontinuum methods. Several quasicontinuum energies have been 
proposed which satisfy the patch test for a limited class of problems, but a quasicontinuum energy 
that satisfies the patch test for many-body potential interactions in several space dimensions with 
general coupling interfaces and coarsening has yet to be developed. 

The quasi-nonlocal energy (QNL) [39J gives an explicit and implemented algorithm for close range 
interactions (up to next-nearest neighbor interactions for a chain and up to fourth-nearest neighbor 
interactions for a FCC lattice) by restoring the symmetry of the interactions. Ghost forces remain 
for nonplanar interfaces and coarsening [19]. The geometric consistency approach gives geometric 
and algebraic conditions to allow finite range interactions, but an implemented algorithm has only 
been given for low index planar interfaces. An efficient and implemented algorithm has yet to be 
given for a general nonplanar interface surrounding a defect. Ghost forces still remain for nonplanar 
interfaces and coarsening [19]. Further, the generalization of the geometric consistency approach to 
three-dimensional problems seems to depend on special decompositions of three-dimensional space 
into congruent tetrahedra (for example, the hexagonal lattice is such a special decomposition of 
two-dimensional space) . 

The QCE energy, on the other hand, can be implemented for any decomposition of three- 
dimensional space into tetrahedra (see (21) in [28])) and so can the blended quasicontinuum energy 
(BQCE) which we propose below. In fact, a code for the BQCE energy can be obtained by simply 
modifying any QCE code and suppressing the coarsening in the blending interface. We have done 
that for our two-dimensional QCE code [40J. The BQCE code will then be able to utilize any 
adaptive atomistic-to-continuum and continuum coarsening features of the QCE code. 

More recently, generalizations have been given for the explicit extension of the QNL energy to 
allow finite range interactions [22,38j. The consistent a/c coupling [38] gives an implemented patch 
test consistent quasicontinuum energy for two-dimensional problems with general pair potential 
interactions, atomistic-to-continuum coupling interface, and coarsening. No generalizations of these 
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or other patch test consistent methods to three-dimensional crystals, multi-lattices, or multi-body 
potentials are yet known. 

Since there remain important problems for which no patch test consistent methods are known, 
we feel it is important to develop general strategies for reducing the error of coupled energies 
which can be applied to a broad class of methods and problems. We propose that by blending the 
atomistic and corresponding Cauchy-Born continuum energies of QCE in an interfacial region with 
thickness of a small number k of blended atoms, a blended quasicontinuum energy (BQCE) can be 
developed with the potential to significantly improve the accuracy of the QCE energy near lattice 
instabilities such as dislocation formation and motion. This blended quasicontinuum energy can 
be implemented for any problem for which a QCE energy has been developed, and we think that 
most existing codes for QCE could be easily modified to implement BQCE. 

The blended quasicontinuum energy is not patch test consistent, but our error estimates in 
Section [3] show that the error due to the ghost force can be reduced significantly by optimizing the 
blending function over a blending region of small size k. Therefore, we think that the use of BQCE 
is a good strategy for reducing the error of QCE when applied to problems for which no patch test 
consistent methods are known. Moreover, patch test consistent quasicontinuum methods tend to 
be more complicated and more difficult to implement than QCE. Thus, even for problems for where 
a patch test consistent quasicontinuum energy is known, we think that the BQCE energy remains 
an attractive quasicontinuum energy. 

The BQCE method is similar to many other methods in which atomistic and continuum energies 
are smoothly blended over an interfacial region, with some features specific to the quasicontinuum 
setting. We call such a region a blending region, and we call such a method a blended method. 
(Other authors call the blending region a "handshake region" , an "interface" , a "bridging region" , 
or an "overlap" .) We call the weights which blend the atomistic and continuum contributions to the 
energy blending functions. For example, the AtC coupling [3], the bridging domain method [6], and 
the Arlequin method for coupling particle and continuum models [5] are all energy-based blended 
methods which are similar to BQCE. By contrast, in other schemes, the atomistic and continuum 
models are coupled abruptly across a sharp interface consisting of only a few atoms. Such schemes 
include the energy-based quasicontinuum (QCE) method [50] . the quasinonlocal quasicontinuum 
(QNL) method [39], the generalized QNL method [22], and the consistent a/c coupling |38j . 

In this paper, we present an error analysis of the blended quasicontinuum energy (BQCE). We 
determine precisely how the error of the BQCE method depends on the size of the blending region 
and how it depends on the blending functions. We analyze the accuracy of BQCE applied to the 
problem of a one-dimensional chain of atoms with periodic boundary conditions and a next nearest 
neighbor pair interaction model. Our choice of a one-dimensional analysis allows us to explicitly 
investigate the accuracy of the BQCE energy for deformations up to lattice instability [16] ; which 
is crucial for the computation of critical strains in lattice statics as well as for the computation of 
the dynamics of defect nucleation and movement. 

We focus on modeling error and do not consider coarse-graining in our analysis. Instead, we 
assume that the continuum energy is a Cauchy-Born energy discretized by piecewise linear finite 
elements with nodes at every atom. We refer the reader to [321 33 i for an analysis of coarse-graining 
for a problem similar to our one-dimensional chain. In addition, we constrain the displacements 
of atoms in the blending region to exactly match the continuum displacement field. In some other 
methods, the displacements of atoms in the blending region are coupled weakly to the continuum 
displacement field using Lagrange multipliers or a penalty method. This is the case in the Arlequin 
method [5], in the bridging domain method [6], and in some implementations of the AtC Coupling 
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0|. Computational experiments comparing various methods of weak coupling are given in ||3.4l6 ll36j . 
A numerical analysis of a mixed finite element method for a weakly coupled problem is given in [5] . 

Blending can also improve the accuracy of quasicontinuum energies which are more accurate than 
QCE, but may not satisfy the patch test for general interactions, nonplanar interfaces, or coarsening. 
We thus also define and analyze a blended generalization of QNL which we call the blended QNL 
(BQNL) energy. The BQNL method is an energy-based, blended coupling which passes the patch 
test for the one-dimensional problem analyzed in this paper, but does not generally satisfy the 
patch test as described above. 

Our analysis extends the techniques developed in [31] for the patch test consistent QNL to obtain 
precise estimates of the error for BQCE, which does not satisfy the patch test. Most importantly, 
our estimates can be used to optimize the blending function and blending interface size of BQCE to 
obtain accurate solutions up to lattice instabilities. Our error estimate in Theorem 15.11 shows that 
the l 2 strain error for the BQCE method can be reduced by a factor of k 3 ^ 2 where k is the number 
of atoms in the blending region. This result is suggested by the results of [I3l[29] on the decay of 
the coupling error for QCE. It has been shown that the QCE method suffers from a 0(1) error in 
the critical strain at which the lattice loses stability (which models fracture or the formation of a 
defect) [16]. We prove in Theorem 14.11 that the error in the critical strain of BQCE can be reduced 
by a factor of k 2 where k is the number of atoms in the blended interface region. In Remark [331 we 
use our modeling error estimates to derive blending functions for the BQCE method which minimize 
the error due to the ghost force. We give our a priori error estimate for BQNL in Theorem 15.21 In 
Theorem 15.21 we show that the t 2 strain error for the BQNL method can be reduced by a factor of 
k 1 / 2 . We prove in Theorem 14.11 that the BQNL method does not suffer from a critical strain error. 

2. An atomistic model and its quasicontinuum approximations 

2.1. The atomistic model problem. We begin by presenting a model of a one-dimensional 
lattice of atoms. Our reference configuration is the integer lattice eZ with lattice spacing e > 0. 
The admissible deformations are N-periodic, mean-zero displacements of uniformly strained states 
of eZ. 

Precisely, for N G N, we let 

JV 

U := {u : Z — > R : = u^+n for all £ G Z and = 0} 

be the set of N-periodic, mean-zero displacements. For F G (0, oo) we let 

4 ■= Fe£ 

be the uniformly deformed state with macroscopic strain F, and we let 

yF '■= {y '■ Z — > E : = + for some u G IA) 

be the set of admissible deformations with fixed macroscopic strain F . We will fix e = -h throughout 
the remainder of the paper, so that the reference length of a period is one, independently of e. For 
notational convenience, we have let the deformations y G be functions of Z instead of eZ even 
though we still think of eZ as the reference configuration. We also remark that the "displacements" 
u which appear in the definition of y& are not displacements measured relative to the reference 
lattice eZ, but are instead displacements relative to the uniformly deformed state y F . 

In the language of mechanics, one would say that our choice of the deformation space 3^f means 
that we have imposed "periodic boundary conditions with macroscopic strain F." Periodic bound- 
ary conditions are mathematically similar to Dirichlet boundary conditions in that constraints are 
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applied to the space of admissible deformations. In our case, the constraint is given by the macro- 
scopic strain F which models the macroscopic scale stretching or contraction of the chain. Periodic 
boundary conditions are useful because they eliminate surface effects near the boundaries of the 
chain. 

For a deformation y G y?, we define a stored energy per period which we call <&(y). Our 
stored energy <&(y) sums the energies of all the interactions between first nearest neighbors and 
second nearest neighbors in the lattice. We compute the energies of interaction using a potential 
4> : (0, oo] — > R which we assume satisfies 

(1) </>GC 4 ((0,oc],M); 

(2) there exists r* so that 4>"{r) > for all r G (0,r*), and 0"(r) < for all r G (r*,oo); 

(3) (j> and its derivatives decay at infinity. 

Examples of commonly used interaction potentials which satisfy the above conditions include the 
Lennard-Jones potential and the Morse potential. Specifically, we define 

JV 

*(y) = ej>( 1 £) + 0(i4 + where ^:= ^~^' 1 . (2.1) 

In (|2.ip . we call the terms of the form 4>{y'^) first nearest neighbor interactions, and we call terms 
of the form <j)(y'^ + second nearest neighbor interactions. 

The reader will observe that we have scaled the energy by e and the interatomic distances by e _1 . 
This is done for two related reasons. First, quasicontinuum methods are only needed when when N 
is large (equivalently, when e is small). By introducing the scaling above, we are able to distinguish 
those parts of the error which are most significant in the practically relevant limit of large N. The 
reader should bear this in mind when interpreting our error estimates in Sections El [H and[5j Second, 
we desire that in the limit as e tends to zero (equivalently, as N tends to infinity), the atomistic 
energy converges to the fully continuum Cauchy-Born energy (|2.5p . This is explained in detail in 
Section \2. 21 below. Roughly speaking, this ensures that the atomistic energy is compatible with the 
continuum energy to which it is coupled in the BQCE and BQNL schemes. The compatibility of 
the two energies is reflected in the fact that the modeling error of the Cauchy-Born energy (|2.4p is 
of order e 2 (See Remark |3.2|) . 

Now suppose we want to model the response of our lattice to a dead load. We let IA* denote the 
algebraic dual of IA, and we equip the space IA with the £ 2 inner product given by 

JV 

< u, v >:= e u^V£. (2.2) 

Generally speaking, a dead load / should be thought of as an element of IA* . However, since U is a 
Hilbert space with the inner product (12.2p . any g G IA* is represented by some /eW. That is, for 
any g G IA* there exists / €Wso that for all v G U, g(v) =< f,v >. Thus, we will think of dead 
loads as elements of U. 

The total energy for an atomic chain subject to a dead load / G IA is then given by 

& otal (y) = $(y)-<f,y>, 
and the equilibrium deformation of the atoms solves the minimization problem: 



y G argmin<3? 

yey F 



(2.3) 
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Let <5$ denote the first variation of the energy $, and let 5 2 Q denote the second variation. We will 
say that a solution y £ yp of the minimization problem (|2.3p is strongly stable if 5 2 <&(y)[u, u] > 
for all u £ U \ {0}, and we will call a deformation y £ 3V an equilibrium if it solves 

5$(y)[u] =< f,y > for all u € U. 

2.2. The Cauchy-Born approximation. We call an energy ^(y) local if it can be written in the 
form 

N 

where ip : (0, oo] — > M is called a strain energy density. Otherwise, we say that an energy is nonlocal 
Observe that energy (|2.ip is nonlocal. The development and use of a strain energy density eft allows 
coarse-graining in the continuum region by utilizing the full range of algorithms, codes, and analysis 
developed for the continuum finite element method. (An alternative quasicontinuum approach 
uses quadrature to approximate the nonlocal coarse-grained energy without the construction of a 
strain energy density \21\ . Peridynamics offer another promising fully nonlocal approach to coarse- 
graining [37].) Thus, it is useful to devise local energies which approximate (|2.ip . The key to 
developing such energies is to observe that as long as y'^ varies slowly between neighboring lattice 
points, we have 

If we replace the second nearest neighbor terms (j)(y^ + y'^+i) by ^{<p(2y'^) + 4>{2y'^ +1 )} in (|2.ip . we 
then obtain the Cauchy-Born energy 

N 1 

^\y) := e£>($ + - {0(2j£) + 0(2|£ +1 )}, (2.4) 

which is a local energy commonly used to approximate (|2.ip . 

An alternative derivation of the Cauchy-Born energy is as follows. Suppose that Y : [0, 1] — > K 
is a C°° deformation of [0, 1], and define y| = Y(e£). It can be shown that 

lim $(y e ) = / V C& ( ^r-(x)] dx =: ^ cb (Y), 
e-^Q+ y [0>1] \dx J 

where ij) ch (r) := (j)(r) + 4>(2r) is called the Cauchy-Born strain energy density. We call the energy 
\E' cb the fully continuum Cauchy-Born energy. We refer the reader to [8j for a detailed discussion 
of the convergence of &(y e ) to ^ cb (Y) as e goes to zero. 

Now suppose that Y is in the Lagrange P 1 finite element space on [0, 1] with nodes at the 
reference position of each atom; so the nodes are e£ for £ = 1, . . . , N. In that case, it is easy to 
show that 

V cb (Y) = $ C V). (2.5) 

Thus, we see that the Cauchy-Born energy <3? cb is a finite element approximation of the fully 
continuum Cauchy-Born energy ty cb . Consequently, we will refer to $ cb as a "continuum" energy. 
We say that the energy is not coarse-grained since & cb has the same number of degrees of 
freedom as the atomistic energy both depend on the deformed position of each atom. 
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2.3. Quasicontinuum approximations. Suppose we want to approximate a solution to the min- 
imization problem (|2.3j) which has defects in localized regions but which is smooth throughout most 
of the lattice. Efficiency requires that most of the lattice be modeled using the Cauchy-Born en- 
ergy (|2.4p . while accuracy requires that the defects be treated using the nonlocal energy (|2.ip . 
Thus, we desire a coupling of the two models. 

We distinguish two approaches to deriving coupled energies. We call these the strain- energy 
based approach and bond based approach. In the strain-energy based approach, one first defines an 
atomistic energy per atom. For our one-dimensional chain, an appropriate energy per atom is 



We call <fr|(y) the atomistic energy at atom £. Observe that is half the total energy of all 

bonds involving atom £, so <3?(y) = e]CgLi ^iv)- The energy per atom should be interpreted as 
the energy per length e at atom £; thus, $^(y) is analogous to a continuum strain-energy density. 

To define a coupled energy we now blend the atomistic energy per atom with the continuum 
energy. Schematically, we choose some blending functions a : [0,1] — > [0,1] and j3 : eZ — > [0,1]. 
Then for a deformation Y : [0, 1] —¥ R with corresponding atomistic deformation y| := Y(e£), we 



In equation (j2.7j) . one should imagine that the function a is zero in a small region containing 
any defects, and one throughout the bulk of the lattice. One should imagine that f3 is one near 
defects, and zero throughout the bulk of the lattice. For an abruptly coupled energy, one should 
imagine that a is the characteristic function of the continuum region, and (3 is the characteristic 
function of the atomistic region. By contrast, for a blended energy, one should imagine that a and 
f3 are more smooth. The energy-based quasicontinuum (QCE) method [28], the bridging domain 
method [6], and the Arlequin method [5] all take essentially the form (I2.7p . Moreover, although the 
AtC method [3] was derived in variational form, the equilibrium condition is the Euler-Lagrange 
equation of an energy similar to (|2.7p . 

We will now explain the QCE energy in detail. Let [0,1] = A U C be a partition of the domain 
[0, 1] into an atomistic region A and a continuum region C. Suppose that Y is in the P 1 Lagrange 
finite element space with nodes at every lattice site, and let y^ := Y(e^) be the corresponding 
deformation of the atomistic lattice eZ. In this special case, the QCE energy reduces to 



(2.6) 



define 




(2.7) 



& ce (y) : = 6 $I(V) + * E where 



(2.8) 



We call the continuum energy at £. 

Our blended energy-based quasicontinuum energy (BQCE) is based on a blend of the atomistic 
energy at £ (|2.6p and the continuum energy at £ (|2.8p . Let 7 : Z — > [0, 1] be a blending function. 
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For y € we define 

N 

^ ce (y) : = 6^ 7 ^|(y) + (1 " 7*)*?(v) 

= c E W + jWM) + Whi)} + We-i + vd + M+i + »£+a)}- 

(2.9) 

We observe that the QCE energy is simply the BQCE energy whose blending function 7 is the 
characteristic function of the continuum region. We will show in Theorem 13.11 that BQCE does 
not pass the patch test where we recall that a quasicontinuum energy passes the patch test if there 
are no forces under uniform strain; for our one-dimensional problem an energy q> cou P led passes that 
patch test if 5<£> C0U P led (y F ) = for all F € (0,oo). We will call a method which passes the patch 
test patch test consistent. 

We will now discuss the QNL and BQNL energies in detail. The atomistic and continuum 
energies of the bond between the nearest neighbors at reference positions e£ and e(£ — 1) are taken 
to be 0(y£). This is the same as the energy of that bond in the fully atomistic model (|2.ip . The 
atomistic energy of the bond between the second nearest neighbors at reference positions e(£ + 1) 
and e(£ — 1) is taken to be (j)(y'^ +1 + y^), and the continuum energy of that bond is taken to be 
\{((>{y'^ + i) + 4>(y'^)}. As discussed in Section [272]. this choice of continuum bond energy is related 
to the Cauchy-Born rule. The coupled energy is then given by 

- 1 
:= 6 £ + e £ + $ + - + (2.10) 

where [0, 1] = A U C is a partition of [0, 1]. 

Our blended quasinonlocal QC (BQNL) energy is based on a smooth blending of the atomistic 
and continuum bond energies used to define the QNL energy (I2.10p . Let 77 : Z — » [0, 1] be a blending 
function. We define 

:= c E W + i&Ue + w^+i) + —^MM) + (2.11) 

We remark that the QNL energy is the BQNL energy whose blending function is the characteristic 
function of the atomistic region A. 

We will now give a simple proof that BQNL passes the patch test. We discuss the modeling error 
of BQNL in more detail in Section [3j The patch test consistency of BQNL is a consequence of the 
following result. 

Lemma 2.1. The set of BQNL energies is the affine hull of the set of QNL energies. In particu- 
lar, any BQNL energy may be expressed as an affine combination of QNL energies with different 
atomistic and continuum regions. 

Proof. The reader may verify that the set of BQNL energies is an affine space. Thus, to prove the 
Lemma it suffices to express every BQNL energy as an affine linear combination of QNL energies 
with different atomistic and continuum regions. Let <$>f l be the BQNL energy with blending 

function j3, and let <&£f be the QNL energy with atomistic region A = {i} and continuum region 
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C = {1, . . . , N} \ {i} for i = 1, . . . , N. We compute 

*g{v) - ® c \y) = e f ^(^ + y'i+i) - hwd) + <f>(2y' i+ i)} 



2 

and so we derive 

N 



*f(y) - ^\y) = e J> Uty + ^ +1 ) - ±{#2$ + *(V t+ i)} 



i=l i=l i=l 



Then we have 



N N 

*fti) = £ Asffo) + (i - E 
i=i i=i 

This expresses ^ q ^ 1 as an affine combination of QNL energies, since we recall that the Cauchy-Born 
energy, $ c6 , is the QNL energy whose continuum region is the entire domain. □ 

In light Lemma 12. 1[ one expects that BQNL will inherit many of the properties of the QNL 
energy. In particular, since any BQNL energy is an affine combination of QNL energies, and since 
the QNL energy is patch test consistent [39J, the BQNL energy passes the patch test. Moreover, 
we show in Remark 14.21 that the BQNL energy predicts the critical strain of the atomistic energy 
as accurately as the QNL method. 

Remark 2.1. (Construction of Patch Test Consistent Blended Methods). Lemma 12.11 suggests a 
general method for constructing patch test consistent, blended methods from patch test consistent 
methods with a sharp interface: one can define a patch test consistent, blended method by taking 
a convex combination of patch test consistent energies with different atomistic and continuum 
regions. Of course, it may be that the method so constructed is not practical. Nevertheless, we feel 
that this observation could be useful in deriving patch test consistent, blended couplings. 

For our analysis, it will be convenient to define a single energy which incorporates both BQCE 
and BQNL as special cases. We will use the blended quasicontinuum (BQC) energy 

N 

*aj(v) ■= e £ <t>(y's) + a^(2y^ + + y^ +1 ) (2.12) 

where a, fi : Z — > [0, 1] are blending functions. We observe that the BQCE energy with blending 
function 7 is the same as the BQC energy with blending functions 

«C : = 7£ == 2 and := 1 - (2.13) 

The BQNL energy with blending function rj is the BQC energy with blending functions 

a? := l _ rj£ := 1 - ^ + ??S ~ 1 and /% := 7£. (2.14) 
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2.4. Summary of notation and auxiliary theorems. Here we collect all the notation and 
auxiliary theorems which we will use below. The reader may feel free to skim this section at first 
and return only as necessary. 
For differences, we write 



Ay ? : = 2/5 - A y f : = y f+1 - 2y ? + 



For means, we write 



2/5+1 - + %£-i - yg-2 



e 3 



._y L +y i _i 
y 5 - 2 ' 



For y G 3V > we will think of y , y , y , and y as N-periodic functions defined on the reference 
configuration. 

We will also use certain bounds on the potential function <p and its derivatives: 

Q(r ) := sup|0«(r)| for i = 1, . . . , 4. 

r>7"o 

We will denote the convex hull of the set A by 

conv A. 

We define the following norms 

{N p 
e £l%l P j for [l,oo), Hi/Hi. : = ^max^|y $ |, 

IMbi.p : = for p e [1,00]. 

Correspondingly, we let l P denote the space U equipped with the norm || • ||^>, and we let U 1,p 
denote the space U equipped with the norm || ■ ||^i, P . Additionally, let Y* denote the topological 
dual of the Banach space Y, and let U~ 1 ' p := (U 1,q )* where p,q G [1, 00] with - + - = 1. We let 
A denote the atomistic region, C denote the continuum region, and I denote the interface. We let 
II ' \\f%(v) De * ne II ' ll^f norm taken over the set V for V = A,C, or X. We denote the closed ball of 
radius r at x in X by 

B x (x,r) := {y G X : \\y - x\\x < r} 

for X one of the spaces £e or IA 1,P . 

We will write 5*ff for the first variation of a differentiable energy functional ^. The first variation 
is a map from yp into U*; we will let 5^(y)[u] denote the first variation of ^ at the deformation 
y G yp evaluated on the test function u G U. We will use the letter u to denote a test function 
belonging to IA throughout the remainder of the paper. The letter y will be used to denote a 
deformation belonging to yp. We warn the reader that in expressions such as 5^(y)[u], u denotes 
an arbitrary test function, not the displacement corresponding to the deformation y. Similarly, we 
will let 5 2 ^ denote the second variation of ^ . The second variation can be interpreted either as a 
map from yp into the space of bilinear forms on U, or as a map from yp into L(JA,U*). We will let 
5 2 ^{y)[u, v] denote the second variation of ^ at y G yp evaluated on the test functions u,v 

We will need the following version of the Inverse Function Theorem which appears as Lemma 1 
in EU. 
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Theorem 2.1. (Inverse Function Theorem) Let X and Y be Banach spaces, let A be an open subset 
of X, and let J- : A — )• Y be a C 1 function. Let xq £ A and suppose: 

(1) \\H*o)\\y <V, and (^(xo)r 1 rnr <*\ 



(2) B x (x ,2 W ) C A; 

(3) || <5.7-"(rci ) — §J 7 (x2)\\l(x,y) — ^ f or °" Xl ' X2 e ^ w ^ H Xl ~~ x 2||x < 27707 

(4) 2Lo 2 n < 1. 

Then there exists x £ X so that J~(x) = and \\x — xq\\x < 2na. 

In Section 5, we will use Theorem 12.11 to prove a priori existence with error estimates for the 
BQCE and BQNL energies. We will apply Theorem 12.11 using bounds on 77 derived from the 
modeling estimates in Section [3] and bounds on a derived from the stability estimates given in 
Section 4. 



3. Modeling Error for the BQC method 

In Theorem I3.lt we estimate the IA~ 1,P norms of the modeling errors of the BQCE and BQNL 
energies. 

Theorem 3.1 (Modeling error in t{ ' p ). Let y £ yp with min^ggy^ > 0. 

(1) Let $!y Ce be the BQCE energy with blending function 7. Recall from equation (|2.13p that 
&j Ce is the BQC energy with blending functions := 7 y£, and /?£ := 1 — 2£±I+2!£lI. y/ e have 

\\6$(y) -6$\ ce (y)\\ u -i, P < Ci||A 2 a € ||, ? + eC 2 ||A/%^||, ? 

+ 6 2 {c 2 ||(i - /%-i)yf|| lf + c- 3 ||(i - my'D 2 \k} 

where Ci := Ci(r mm ), i = 1,2,3, for r mm := 2 min^^y^- 

(2) Let <I>^ n ' be the BQNL energy with blending function /3. We have 

\\6$(y) - 5<S>f\y)\\ u -^ < eC 2 ||A/%^||, ? 

+ e 2 {C 2 ||(l-/3 5 _ 1 )yf||, ? + C'3||(l-/3 e )(^) 2 lk} 
where Q := Ci{r min ), i = 2, 3, for r min := 2 min ?eZ y' c 



(3.1) 



(3.2) 



Proof. We begin by writing down the first variations of <& a g and <£. For y G yp and u € U 1,v , we 
have 



iV 



8^{y)[u] = e <t>'{y'i)u'z + 2a e 0'(2^)4 + /3^'(y^ + y' H1 )(u^ + u' c+1 ) 
5=1 

N 

5=1 

The first variation of $ is a special case of the above. We have 

N 

8${y)[u] = e YWiy't) + + $ + W« + ^+i)K" 

5=1 
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Using these formulas we compute the modeling error 
M[u] : = 6$(y)[u] - 8$ a ,p(y)[ U ] (3.3) 

N 

= e [{(i - Pe-iWt-i + vd - + {(i - m'iy't + ^+1) - 4 

Now we expand the terms 

G e := {(1 - /3 e _ 1 )^'(^_ 1 + $ - c^'(2y£)} and ^ := {(1 - /%)<^ + - a^(2^)} 
in (|3,3p at 2y'^ using Taylor's theorem. We obtain 

e_ 2 
2 



G c = (1 - " « c )0'(2^) + (1 - /3^ 1 )(-60"(2^)^_ 1 + ^(Oi,^.!) 2 ) and (3.4) 



^ = (1 - /% - a*) W{) + (1 " PtW(MW{ + f 0'"(O 2 ^)(y^) 2 ), (3.5) 

where Oi^ E conv{2y^, + y^}, and C>2,£ G conv{2y^, + y^}. When we substitute (13.4|) 
and (|3.5p into (|3.3p . we obtain 



iv 

MM = e ^ 2{1 - a e - ft^WeK 

+ e {(i - PzWWtWt - (i - &-i)^W*)i£-iK (3 ' 6) 



+ £{(i - fo^lPuM? + (i - ft-i)^(Oi^)(i^-i) 2 K- 



Expanding the term of order e on the right hand side of (|3.6p . we compute 

AT 

M[u\ = e 2{1 - « C - ^}0'(2^)4 

+ {-eA/3 e </>"(2^' + e 2 (l - ft-i^W*)^" }u' 6 (3 " 7) 



+ - ^)r(o 2 ,0(y'i) 2 + (i - A^)(»?-i) a K- 



We recall from equation (|2,14p that the BQNL energy f&'g 11 with blending function j3 is the BQC 

energy <3?^ with a>£ = 1 — In that case, we compute that for the BQNL energy <3?^ n ', the term 
of order zero in equation f|3.7j) vanishes, and so 

N 

+^{(1 - /%)<(0 2 , e )(y£) 2 + (1 - ^-i)r(Oi,e)(^-i) 2 K- 

Therefore, by Holder's inequality, 

|H>(y)[u] - ^(y)M| < {eC 2 ||A/3^||, f + e 2 C 2 \\(l - ^-M'\k + ^11(1 ~ ft)(^) a ||«} Kll* 

where o := -^r. Thus, 
^ p— l ' 

||<5<%)[u] - 6*f (y)[u]\\ u -i, P < eC a ||Aft»?||* + e 2 {C 2 \\(l - ||, ? + C 3 ||(l - /%)(^) 2 |k}- 
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This proves the first claim made in the statement of the theorem. 

On the other hand, for the BQCE energy <3?^ ce there is a ghost force arising from the term of 
order zero on the right hand side of equation (I3.7p . Using formulas f)2. 13f) . we compute 

2{1 — a(. — /%} = a£ + i — 2a% + ct£-i = A 2 a^. (3.8) 

Substituting this expression in equation (13.7p . we derive 

N 

5<S>(y)[u] - 5^ ce (y)[u] = e J> 2 a^W € K 

+ {-eAWW^ + e 2 (l - Pt-iWWtWH 

+ ^{(1 - ^)<t>"\O u ){yl? + (1 " Pt-i)<t>"'(Oi,s)(yU) 2 H- 

Thus, 

- <^ ce (y)[u]| < /ciHA 2 ^^ + eC 2 ||A/3^'||, f 

+ e 2 C 2 \\(l - Ps-i)y?\\ e p + e 2 C 3 \\(l - Ps)(yl)%p\\\u'\y e . 

This proves the first claim made in the statement of the theorem. □ 

Remark 3.1 (BQNL is patch test consistent and BQCE is not patch test consistent). Estimate fj3.2[) 
implies that BQNL is patch test consistent. For observe that if y F is the uniform deformation, then 
(y F )'l = (y F )'l' = for all £ £ Z. Thus, by ||^(^)-^" i (^)|| w - 1 . P = ^ V)lb^ = 0- 

On the other hand, observe that the term Ci||A 2 ag||^p which appears on the right hand side of 
estimate (|3.ip does not vanish under uniform strain. This reflects the fact that BQCE is not patch 
test consistent. 

Remark 3.2 (Interpretation of modeling estimates). We will now give a detailed interpretation of 
the each term in estimates (13. 2D and (I3.ip . First, we consider the term 

e 2 {C 2 \\(i - Pz-^'IUp + a 3 ||(i - frWifk}- 

Observe that the function 1-/3^ is supported in the interface and continuum, and that it is 
identically one throughout the continuum. This suggests that the term arises from the error of 
the Cauchy-Born model. Recall that the Cauchy-Born energy is the BQNL energy with blending 
function /3g = 0. Therefore, by the modeling estimate for BQNL given in Theorem 13. II we have 

\m V )-6#*(3,)\\u-i* < ^{C2\\y'l'\k + C 3 \\(y'l)\ P( }. 

Consequently, we will call the term e 2 {C l 2||(l — /?^_i)y^"||^p + C* 3 1 1 ( 1 — /%)(y£') 2 ||^p} the Cauchy-Born 
error. 

Next, we consider the term eC^HA/^y^H^. We observe that A/3^ is supported in the interface. 
Thus, the term eC^HA^y^H^p arises from the error caused by coupling the atomistic and continuum 
models in the interface. We will call this term the coupling error. 

Finally, we consider the term Ci || A 2 a^||^p which appears in the modeling estimate for BQCE but 
not in the estimate for BQNL. Let y F be the uniform deformation y F := Fet;. We call 5^ ce (y F ) the 
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Continuum 
* IntertacG 
o Atomistic 



Figure 1. Graph of the blending function 7 with interface region X = J\ U Jz- 

ghost force associated with the energy Qj, and we observe that under uniform strain formula (I3.7P 
reduces to 

N 

8^ ce (y F )[u] = e^-A 2 a^(2F)u' v 

Thus, we see that <p'(2F)\\A 2 a^\\iP is the IA~ 1,P norm of the ghost force. We will call C\ \\ A 2 a^\\pv 
the ghost force error. 

Remark 3.3 (Dependence of ghost force error on interface size). We will now analyze the dependence 
of the ghost force error on the size of the interface and the shape of the blending function. Our 
analysis will lead to an estimate for the rate at which the ghost force error decreases with the 
number of atoms in the interface, and to an optimal family of blending functions for the BQCE 
method. Consider the BQCE energy $^ ce with blending function 7 as depicted in Figure [TJ 

Each transition between the atomistic and continuum models results in some ghost force. We 
will consider the ghost force which arises from a single transition from the atomistic model to the 
continuum model. The total ghost force is, of course, the sum of the ghost forces due to each 
transition. Let J2 ■= {i, ■ ■ ■ , i + k} be part of the interface, which we will denote simply as J is the 
following. Assume that atoms i — 2, i — 1, and i are in the atomistic region, and that atoms i + k 
and i + k + 1 are in the continuum. That is, suppose that 7(£) = for £ € {i — 2, . . . ,i}, and that 
7(£) = 1 for £ G {i + k, i + k + 1}. We will estimate the ghost force due to the transition which 
occurs over region J . 

First, we construct a blending function which implements such a transition. Let 7 S : [0, 1] — » [0, 1] 
be a twice continuously differentiable function such that 7 S (0) = 0, 7 S (1) = 1, and 7^(0) = 7^(1) = 0. 
Extend 7 S to a function defined on R by taking J s (x) = for x € (— 00, 0) and 7 s (x) = 1 for 

x G (1, 00). Now let J+ := {i - 2, . . . , i + k + 1}, and define jj : J+ -> [0, 1] by 7t7 (f ) = 7s . 

We call 7 S the shape of the blending function 7^7. 
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We will show 



2 - i— 2 

l A lj\\l P {J) < epkp 



dx 2 



LP ([0,1]) 



for all p G [1, oo]. Suppose p G [1, oo), and compute 



ll A2 7j||ff(j) = e 2 \Wj\\£P(j) =e p kp 2 < 
By Taylor's theorem, 



1 fc 
f=0 



7s 



^)-27 s (|)+7 s (^ 



l 



(3.9) 



7s i k 



27s (! 



+ 7s(V 



i-i cix 2 



(s)rj 6 (s)ds 



(3.10) 



where 



%0) :r 







+ if s € 

otherwise. 



fc ' fc 



We note that (|3.1U|) holds for £ = i and t; = i + k only since we have assumed that 7^(0) = 7^(1) = 0. 
Now observe that rj^ is a non- negative function whose mass is one, and that x *— > \x\ p is convex for 
p > 1. Therefore, Jensen's inequality implies 



< 



<ix 2 



00 



This yields the estimate 



I a2 7^II^(^) - epA;p 



epkp 



-2 



dx 2 

p k 



(is 



<ix 2 



epkp 



-2 



As 



dx 2 



LP ([0,1]) 



for all p 6 [1, oo). By a similar argument, one can show 



dx 2 



A2 lj\\e°°(j) < k 2 

"i°°([0,l]) 

Recall from Remark 13.21 that the ghost force error due to the transition over J is 

2„, II _ ._ II A 2- 



ll A aj\\eP(j) ■= ||A 7j _ ll^f( > 7) 
By Minkowski's inequality and the estimates above, we have 



--2 



d 2 ls 



dx 2 



(3.11) 



l A a Aez{j) = ll A uuuj) ^ ll A 7^ll^(j) < t p kp~ 

for all p £ [1, oo]. This gives the dependence of the ghost force error on e, the blending function, and 
k. We remark that inequality (|3.11j) is sharp. In fact, if j s is sufficiently smooth, then HA^yH^vj-) 



converges to epkp 



dx 2 



LP ([0,1]) 



as k tends to infinity. 
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We pause here to explain why we assumed 7^(0) = 7 S (1) = 0. It can be shown that no estimate 

- -—2 

of order e^k? holds unless 7 S (0) = 7 S (1) = 0. We leave the proof of this fact to the reader, and 
will instead give an illustrative example. Suppose that j s were the linear polynomial with 7s (0) = 
and 7 S (1) = 1, so 7 S is not differentiable at or 1. Then 

= zhh- 1 . (3.i2) 

In Section [51 we use an estimate of the IA~ 1 ' 2 norm of the modeling error in order to obtain 
convergence results. Thus, we are particularly interested in the U~ 1,2 norm of the ghost force error. 
By estimate (|3.1ip . if 7 S satisfies 7 g (0) = 7 S (1) = then the U~ 1,2 norm of the ghost force error 
decreases with k as k~2 . On the other hand, if 7 S is the linear blending function then we see that 
the ghost force error decreases as k~ . We conclude that if a large blending region is desired, it 
is best to choose a blending function which satisfies j' s (0) = j' s (l) = so that the faster rate of 
decrease is obtained. In particular, we suspect that the linear polynomial 7 S is a poor choice for 
the shape of the BQCE blending function. 

We will now use estimate f)3.11j) to derive an optimal shape for the blending function of the 
BQCE energy. As discussed above, we will use the U~ 1,2 modeling estimate to prove our error 
results in Section Thus, we would like to find a family of blending functions which minimizes the 
U~ 1,2 ghost force. Estimate (|3.1ip suggests that we should choose the shape 7 S which minimizes 



subject to 7s(0) = 0, 7s(l) = 1, and 7 S (0) = 7 S (1) = 0. The Euler-Lagrange equation 

L 2 ([0,1]) 



dx' 1 

for this minimization problem is 



(x) = for all x € (0, 1) 



dx 4 

Thus, the optimal shape 7 S is the cubic polynomial which satisfies the constraints 7 S (0) = 0, 
7s(l) = l, and 7^(0) =7^(1) = 0. 

Remark 3.4 (Dependence of coupling error on interface size). Now we examine the dependence of 
the coupling error on the number of atoms in the interface and the shape of the blending function. 
First, we consider the coupling error of the BQNL energy. Following the notation established 
in Remark 13.31 let J := + k} be contained in the interface. Let j3 s : [0, 1] — > [0, 1] be 

a continuously differentiable function with /3 S (0) = and (3 s (l) = 1. Extend j3 s to a function 
defined on M by taking j3 s {x) = for x E (— oo,0) and /3 s (x) = 1 for x 6 (l,oo). Now let 
(3 j : {i, . . . ,i + k} — > [0, 1] be defined by = /3 S ( j?)- Using an argument similar to the one 

given in Remark 13.31 one can show 

WA^^c^k^- 1 

for all p G [1, oo]. The result above yields the estimate 

z\\^Pjy"\\e e {j) < 4 A Pj\\e p c (j)\\y"\\e^(j) < ^^k^c^yw^^y 

This gives the dependence of the coupling error of the BQNL energy on e, the blending function, 
and the size of the interface. 

We now address the coupling error of the BQCE method. Let the blending function, jj, of the 
BQCE energy be defined as in Remark 13.31 We recall from equation fj2. 13f> that the BQCE energy 
with blending function "fj is the BQC energy, ^ a ,/3, with 
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Thus, using Minkowski's inequality and the result for the BQNL energy, we make the estimate 



4^Pjy"\\e t (j) ■-- 



A(1 _ 7g + l+7g- 1 )/ 



< 4 A ljy"\\e?(j) < e l+p kv 1 C 7 ||y"||^ c( i7 ). 



(3.13) 

This gives the dependence of the coupling error of the BQCE method on e, the blending function, 
and the size of the interface. 

Remark 3.5 (Higher order estimates). Suppose we let the measure of the interface be a fixed fraction 
v of the total measure of the domain as e tends to zero. Then the number of atoms k in the interface 
would be approximately ^, and estimates fj3. 11 f) and (I3.13P would reduce to 

II a2 «^I!^U) < Ce 2 uv- 2 and e|| A^jy"]^^ < De 2 vv~ l . 

Remark 3.6 (Newton's Third Law). We observe that the forces arising on each atom due to the 
BQC energy can be decomposed into a sum of central forces which satisfy Newton's Third Law. 
In particular, the forces arising from both the BQCE and BQNL energies satisfy Newton's Third 
Law. 



(4.1) 



4. Stability of the BQC method 

First we derive expressions for the second variations of the atomistic and BQC energies. We 
have 

N 

5 2 ^(y)[u,u} =eY,r(y'M\ 2 + 4a^"(2^)|^j 2 + ptf'ty + y^ +1 )\u^ + u^ +1 | 2 

JV 

= eZ>V*)Kl 2 + 4c^"(2^)|4j 2 + &0Ve + ^+i)[ 2 Kl 2 + 214+11 - ' "t 
5=1 

N 



2 ,2|,//|2i 

— c lit 



where 



{y'i + y'z+i) + + yd + «^"(2^) 



and 



(4.2) 



2' ST w? " s_rJ -' 2 

S € :=-^Wf + l4+i)- 
The second variation of the atomistic energy is, of course, a special case of the above. We have 

* ri ii 

5 2 d>(y)[u, u] = e £ 0"(^) + 4 + ^ +1 ) + + y^)) \u^\ 2 + (-eW* + ^ +1 )|4f) 

5=1 1 J 

AT 



e ^^|4| 2 + e 2 ^|4'| 
5=1 



where 



As := 0"(^)+4 



and 



(4.3) 
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We begin our analysis with a lemma bounding \Ap — A^\. In Remark 13.21 we interpreted each 
term which appears in the estimate below, and in Remarks 13.31 and 13.41 we explained how each term 
depends on e, the blending function, and the number of atoms in the interface. We concluded that 

ellA^T/gl^oo < ek^y'W^^ and \\^ 2 a^\\^ < k~ 2 , 

where k is the number of atoms in the interface. The reader should keep these scalings in mind 
throughout Section [H 

Lemma 4.1. 

(1) Let <l?^ ce be the BQCE energy with blending function 7. We have 
max |yig — A^\ <2C , 2||A 2 q^||^oo 

+2eC-3||A%£'||*» + 2e 2 {C 3 ||(l - ft_i)j/£V + C 4 ||(l - &)(l/£) 2 |l*»} 

where a ? := 7$, := 1 - 22±1+Iti j and q = d{r min ), i = 2, 3, 4, for r min := 2 min ?eZ y£. 

(2) Zei $^ n ' be the BQNL energy with blending function /3. We have 

max|^ - ^| < 2 e C 3 [|A/^'||^ + 2e 2 {C 3 ||(l - + C A \\(1 - /%)(y£) 2 ||l~) 

where C{ = Ci{r mm ), i = 3,4, for r mm := 2 min^ e zy^. 

Proof. The proof of the lemma is extremely similar to the proof of our modeling estimate above. 
For the general BQC energy $ Qj/ g we compute 

A^-A^ = 2(# - l)0"(y£ + y' t+l ) + 2(/%_i - + y\) + ^4>"{2y\) 

using formulas (|4.2p and (|4.3p . Then we expand all terms above at 2y^ using Taylor's theorem. So, 

Af-At = 4{^ + a 5 - 1}^'(2^) 

+ 2e<T(2^){(/% - l)y£ - - 1)^} (4.4) 

+ 6 2 {(/3 5 - l)0( 4 )(0 2 , c )(y£) 2 + (fc-i " IJ^CO^Jd^-i) 8 } 

where Oi^ S conv{2y'^,y'^_ 1 + y^}, and 2 ,£ £ conv{2y^, y^ +1 + y^}. Now we expand the second 
term in curly braces on the right hand side of (|4.4p . We obtain 

A\ - A e = 4{Pz + a e - l}0"(2y£) 

+ 2{eA^y£ + e 2 (/3 5 _ 1 - l)yf }^'"(2^) (4.5) 

+ e 2 {(/% - l)«^ 4 )(0 2 , e )(y^) 2 + (ft-i - l)0 (4) (Oi >? )(^-i) 2 }- 

We now consider the case of the BQCE energy. As a consequence of equation (|3.8p . we have 

4{/% + a 5 - 1} = -2A 2 a 5 . (4.6) 

Thus, by substituting (|4.6p into (|4.5p . we see that for the BQCE energy 

|^ - A e | < 2C 2 ||A 2 a 5 ||,oo + 2eC 3 ||A / 9 CI 4V + 2e 2 {C 3 ||(l - /%_i)yf ||^> + <7 4 ||(1 - /%)(y£ ) 2 ||^}. 

This proves the first claim made in the statement of the lemma. 

On the other hand, we recall from the proof of Theorem 13.11 that for the BQNL energy, 

/% + a € - 1 = for all £ G Z. 
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Therefore, for the BQNL energy, the first term in curly braces in equation (|4,5p vanishes, and we 
have 

|^ - ^| < 2eC 3 ||A/3^|| £ oc + 2e 2 {C 3 ||(l - ft-i)j£"||*» + C 4 ||(l - /%)(^') 2 |l^}- 
This proves the second claim made in the statement of the lemma. 

□ 

Now we derive estimates relating the U ,2 coercivity constants of the Hessians of & a ,/3 and <E>. 
Define 

c(y)= inf 6 2 ^y)[u,u], cp(y) = inf 8 2 ^ nl (y)[u, u], 

\\ u ' »2=1 \W ,2=1 

c 7 (y)= inf 5 2 ^ ce (y)[u,u}. 

u' ,2=1 

For our a priori error estimate, we would like to show that if y a is a strongly stable minimizer of 
<&, then cp{y a ) > c(y a ) and c 7 (y a ) > c(y a ). However, such a general result was not proved for the 
QNL method in [31]. Instead, a weaker stability result that is restricted to "elastic states" without 
defects |32] was proved. We now extend this a priori stability result to the BQC method. 

Remark 4.1 (A bound on second neighbor interactions). We will place an additional condition on 
the set of admissible deformations in order to prove our a priori and a posteriori stability estimates 
in Theorems 14.11 and 14.21 We will assume that 



v 

niin^ > —. (4.7) 



Under this assumption, the constants and in the expressions for the atomistic and BQC Hes- 
sians are nonnegative. Assumption (14. 7ft is justified since y'^ < ^- only under extreme compression, 
and in that case the second nearest neighbor pair interaction model (|2.ip itself can be expected 
to be invalid. The authors of [31], [13], and [12] all consider energies similar to ()2. 1 1) . and they all 
make assumptions similar to (|4.7j) (see Section 2.3 of [31J for further discussion of this point). 

Theorem 4.1 (A priori stability). Let A = ming A^(y). Assume that ming y'^ > ^- > 0. 

(1) Let <3? 7 ce be the BQCE energy with blending function 7. We have 

c 7 (y) > A-2C 2 \\A 2 a i \\ioo 

-2eC'3||^^ / ||£ o - 2 e 2 {a 3 ||(l - ^-O^'IU- + (^411(1 - ^)(2yg) 2 H^}, 

where a c := 7^ ft ■-l-^±lp£=l j and Ci = Ci{r min ), i = 2,3,4, for r min := 2 min ?eZ ^. 

(2) Ze£ 6e f/ie BQNL energy with blending function /3. We have 

cp{y) >A- 2eC 3 || Afty£|| € ~ - 2e 2 {(7 3 ||(l - ft-Oyf ||/« + C 4 ||(l - ft)(y£) 2 |l^}, 
where Ci = Ci{r mm ), i = 3,4, /or r mm := 2 min^^y^- 
Proof. Let u G ZY ' 2 with ||u||^i,2 = 1. Recall that by formula ()4. 1[) we have 

N 

5 2 ^(y)[ U ,u} = e^^l^l 2 + e 2 ^|4'| 2 
5=1 
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for any BQC energy ^ a ,/3- Since we assume min^ y'^ > all the coefficients are nonnegative, 
so 

N 

5 2 $ a ^{y)[u,u} > e^i e |4| 2 . 

Now we assume that the energy <& a ,p is a BQCE energy. Then using the first part of Lemma 14.11 
we have 

N 

8 2 ^ ce (y)[u,u] > e^A^l 2 - 2C 2 || A 2 a € || £ oo 

- 2e<5 3 ||A%£||*» - 2e 2 {C 3 ||(l - /%_i)yf ||*» + C 4 ||(l - ^W{) 2 \\tPo} 
>A-2C 2 \\A 2 

- 2e<5 3 ||A/3^'||*» - 2e 2 {C 3 \\(l - ||*» + C 4 ||(l - ^)(^') 2 |l^}- 

This proves the first claim made in the statement of the theorem. 

For the BQNL method, we make a similar estimate using the second part of Lemma 14.11 We 
have 

N 

6 2 $f(y)[u,u} >e^A^\ 2 - 2eC 3 ||A/%/£|k~ - 2e 2 {C 3 ||(l - + 04(1 - /%)(^) 2 |l^} 

> A - 2eC 3 ||A^^||^ - 2e 2 {C 3 ||(l - /VO^'lk- + C 4 \\(l - ^)(y'l) 2 \\^}. (4.8) 
This proves the second claim made in the statement of the theorem. □ 

Remark 4.2 (Accuracy of critical strain). Let y F be the uniform deformation y F := Fe^. Let F* 
be the strain at which the lattice loses stability in the atomistic model, so 

F* := inf{F G (0, oo) : c(y F ) < 0}. 

We call F* the critical strain of the atomistic model, and we define the critical strains of the Cauchy- 
Born, BQCE, and BQNL energies similarly. The significance of the critical strain is discussed in [16]. 
Under uniform strain, the U~ 1,2 coercivity constant c cb (y F ) of the Cauchy-Born energy is 

c ch {y F ) = <P"{F)+^"(2F) =:A F . 

Furthermore, by formula (|4.3p . for the atomistic model under uniform strain A : = min^^ = Ap. 
Thus, by Theorem 14. 1\ we see 

c 7 (y F ) > A - 2CyA 2 a £ |U = c cb (y F ) - 2C 2 ||A 2 a £ ||^o > c cb (y F ) - 2C , 2 C 7 A;~ 2 , and 

(4 9) 

cp(y F )>A = Ccb (y F ), 

where C 7 is the constant which arises in estimate (|3.1ip . and k is the number of atoms in the 



interface. In the language of [16] , estimates (j4.9j) imply that there is no error in the critical strain 
of the BQNL energy, and that the error in the critical strain of the BQCE energy decreases as k~ 2 . 

For an a posteriori existence result, one would like to show that if y bqc is a strongly stable 
minimizer of $ a ,p, then c(y bqc ) > eg(?/6q C ). We will not present an a posteriori existence result, but 
we give an a posteriori stability estimate anyway. Using this estimate, one can prove a posteriori 
existence theorems similar to Theorems 15.21 and 15.11 



Theorem 4.2. (A posteriori stability) Assume that ming y'^ > ^- > 0. 
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(1) Let <l? 7 ce be the BQCE energy with blending function 7. We have 
c(y) > c 7 (y)-2C , 2 ||A 2 a ? ||£oo 

-2e(7 3 ||A/%y£[k~ - 2e 2 {C 3 ||(l - + C 4 ||(l - /%)(y£) 2 ll*~}- 

w/iere a ? := 75, /% := 1 - ; and A = Q(r mm ), i = 2,3,4, /or r min := 2 min ?eZ ^. 

(2) Let §f l be the BQNL energy with blending function ft. We have 

c{y) > c p (y) - 2eC 3 ||A/%y£||^ - 2e 2 {C 3 ||(l - ^_i)yf ||*» + C 4 ||(l - Pt)M[) 2 \\t°°}, 

where C t = d{r mm ), i = 3,4, for r mm := 2 rnin eeZ y£. 

Proof. The proof of the theorem is extremely similar to the proof of Theorem l4.1i We use Lemma [4.1l 
and that B^ > whenever min^ y'^ > ^- . □ 

5. A PRIORI ERROR ESTIMATES FOR THE BQCE AND BQNL METHODS 

We are now ready to prove our a priori error estimates. Our results generalize Theorem 8 and 
its proof given in |31j . Before stating the theorem, it is convenient to establish some notation. 
Let y be a minimizer of the energy <& total . We will use Theorem 14.11 in the proof of Theorem 15.21 
Thus, we must assume that y is an elastic state. That is, we assume that A := ming A^ > for 
5 2 &(y). See the discussion preceding Theorem 14.11 for more explanation. Now define the energy 

Qtotal . y p _j. R by $total( y } ._ $**(y)_ < f,y >. 

We realize that the statement of Theorem 15. II may seem slightly complicated, but the basic idea 
of the theorem is quite simple. In essence, Theorem 15.11 states that if y is a stable minimizer of the 
atomistic energy which is sufficiently smooth in the continuum region, and if the ghost force error 
is sufficiently small, then there exists a stable minimizer y 7 of the BQCE energy which is close to 
y. The conditions f)5. 1 j) and (|5.2p make the hypotheses that y must be sufficiently smooth in the 
continuum region and that the ghost force must be small precise. 

Theorem 5.1 (A priori error estimate for the BQCE method). Let &° tal } y > and A be as above. 
Assume that A > 0, and that min^y^ > ^- > 0. There exist constants 5i := o~i(A) and 5 2 '■= 
^(min^ y'^A, Ci, C 2 , C3, 7) so that if 

2C' 2 ||A 2 a 5 ||^ +2eC 3 ||Aft 1 4V + 2e 2 {C' 3 ||(l - /%_x)yf ||^ +C 4 ||(1 - /%)(y*p 2 |M < <*i, (5-1) 
and 

e-sCiHA 2 ^^ + e^C 2 \\A^y'l\\ ei + e§ {C 2 ||(l - /3^i)yf \\p e + C 3 ||(l - /%)(^') 2 |M < h, (5.2) 
then there exists a locally unique minimizer y 7 of &° tal which satisfies 

||y-y 7 Li, 2 < j [Ci\\A 2 a^y e +eC 2 \\Ap^y e + e 2 {C 2 \\(l -/%_x)yf + C 3 \\(l - (3^)%,}} . 

Proof. Let T : U 1 ' 2 -> U^ 2 by F{w) := 5&° tal {y + w). We will apply Theorem O to T with 
xq := 0. In order to apply Theorem 12.11 we need to find a bound 77 on the residual ^(O), a bound 
a on || ((^(O)) -1 112/^-1,2^1,2), and a Lipschitz constant L for 5F on the ball Byi,2(0,2r]a). We will 
find 77 using the modeling estimates given in Theorem 13-H and we will find a using the stability 
estimates given in Theorem 14.11 together with the condition (|5.ip . Once we have these constants, 
we must verify the condition 2Lo~ 2 rj < 1. As we will see, this condition holds if 6\ and 5 2 are taken 
sufficiently small. Once these facts are established, Theorem 12.11 implies that there exists y 7 G ^f 
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such that J~"(y 7 ) = and \\y — y 7 ||^i,2 < 2r/c Finally, we show that y 7 is a strongly stable minimizer 
of <^ t ° tal which proves Theorem 15.21 

1. Modeling error: ||.F(0)||^-i,2 < rj. First, we find r\ so that ||.F(0)||^-i,2 < rj. By Theo- 
rem [37TJ we have 

||^(0)|| w - ll2 = \\5^(y)\\ u - lt , 

< Ci\\A 2 a^ + eC 2 \\A^y'l\\ Pe + e 2 {C 2 ||(l - ^-i)^^ + C 3 ||(l - ^)(^) 2 |l^}, 

so we take in Theorem 12.11 

r, := Ci||A 2 a^2 + eC 2 \\A(3^\\ £ 2 + e 2 {C 2 ||(l - Ik + ^IIU - hMYWe*}- 

2. Stability: |] (<5^ r (0)) — \\utl- x '' 1 ,u x ' 2 ) — °"- Second, we find <r so that 

II (^(o))- 1 !!^-!,^) <(7. 

By Theorem 14.11 we have 

cp(y) > A-2C 2 ||A 2 q 5 ||^ -2eC 3 ||A/3^||^ - 2e 2 {C 3 ||(l - /3 5 -i)^'||^ + ^([(l - /3 5 )(^) 2 |k-}- 

(5.3) 

Then if we take 6% < =j , we can combine (|5,3p and (|5.ip to find 

c/s(y) > A " 2C 2 1| A 2 a ? ll.oo - 2eC 3 1| A/%^' - 2e 2 {C 3 ||(l - + ^||(1 - /%)(y£) 2 |W 

> A-tfi 
A 

~ 2"' 
In that case, 

1 ,2 
c/jfo) " I' 
So we can use 

a := |. (5.4) 

3. Lipschitz bound of 5 J- on 5^1,2(0, 2rja). Now we bound the Lipschitz constant of 5J- on 
S^i,2(0, 2^cj). The Lipschitz bound follows easily if we can assume that 

^ + <4>~^ (5.5) 

for all u> with ||u;|| w i,2 < 2r/cr. We will choose 5 2 so that (|5.2p implies (|5.5p . To that end, observe 
that for 1 1 w | |^i,2 < 2t/<t we have 

/II —ill /II — i / — A 

||w lk°° < e 2 1|^ 11^ < 2e = M x e 277 

where M[ := 2a. Then we need only choose <5 2 < ^gr min^ y'^ to ensure that ()5.5|) holds. 
To see the Lipschitz bound, let w\,W2 £ B^i, 2 (0, 2r?<r). Then expand 

|£ 2 ^ ce (y + «>!)[«, u] - <5 2 ^ ce (y + w] | 

using formula (|4.ip . One sees easily that 

la^Hy + wOM -£ 2 <M ce (y + ^)M < l'\\ w [ - w ^\\u'\u\\v'\\ g 
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where V := L'{C^{^ min^ y£)). In that case, we have 

\\5F{ Wi )-5Hw2)\\l { u^m-^) = \\5 2 ^ ce (y + w 1 )-5 2 ^ ce (y + w 2 )\\ L{ul ,2, u ^) 

— 1 /n 

< e 2L \\wi — W2\\ui,2 = ^H^i ~~ ^ ll^i.a 

where L := eT^U . 

4. The fixed point condition: 2La 2 n < 1. We show that if #2 is sufficiently small then 
2La 2 n < 1. We have 

_i 1 -4 2 

2La 2 r? < 1 <^ 2er2L'a 2 r) < 1 e~2r/ < (5.6) 

8i 

Observe that (|5.2p could be restated as e _ 5^ < <5 2 with 7/ as defined in part (1) above. Thus, 
inequality (|5.6|) holds whenever 62 < -gjj- 

5. The error estimate. Take S\ := ^, and 62 < min{jj7, 2 w ( }■ Then by Theorem 12.11 and 



the conclusions of parts (1-4) above, there exists y 7 E U 1,2 with 5&° tal (y^) = and 



-7 

Hz/7 - yllw 1 - 2 ^ 2r i a - 

It remains only to show that the equilibrium y 7 is a locally unique minimizer of <& t ° tal . We will 
show that 5 2 ^ ce (yg) is positive definite. We have 

c 7 (y 7 ) > c 7 (y) - |c 7 (y) - c 7 (y 7 )| >~~ L(2rja) > - - - = 0. 
The second inequality follows from (|5.4p . the Lipschitz bound, and the error estimate 

IK - y\\uw < 2 w- 

The third inequality follows from the condition 2La 2 r] < 1. □ 

We will now restate Theorem [5J] using the estimates given in Remarks l3.3l and l3.41 For simplicity, 
assume that all transitions between the atomistic and continuum models occur over regions which 
contain k atoms. Let a and (3 be the blending functions so that $ 7 ce = $ a g. (See (|2.13j) for the 
precise definitions of a and /3 given 7.) Using the estimates given in Remarks 13.31 and 13.4} let Cj 
be a constant so that 



e||A/%|| €2 < Cjeh-^ and e||A/% \\ eoo < C{ek-\ 



and let C 2 be a constant so that 

||A 2 a ? || £ 2 < ClehhT^ and ||A 2 a $ ||^ < C^AT 2 . 

In essence, Corollary 15.11 states that if y is a stable minimizer of the atomistic energy which is 
sufficiently smooth in the continuum region and if the ghost force error is small, then there exists 
a minimizer y 7 of the BQCE energy (^ t ° tal which is close to y. 

Corollary 5.1. Let <& l ° tal } y ) and A be as in Theorem \ 5.1\ and let C\, Cr[, and k be as above. 
There exist constants Si and 62 so that if 

2k- 2 ClC 2 + 2ek- l ClC ? \\yl\\^ {x) + 2e 2 {C 3 ||(l - ||*» + C A \\{\ - /%)(^') 2 |K} < <*i, 

and 

hT*CqC x + efc-ic7(7 2 ||^ / ||, 00(x) + ei {C 2 ||(l - P^y? ||, ? + C 3 ||(l - PtWtfk) ^ 5 ^ 
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then there exists a locally unique minimizer y 7 of &° tal which satisfies 



\\y - y-rWu 1 ' 2 

4 

< — 
- A 



eh-l (%& + el k-^aiC 2 \\y'l\\^ {x) + e 2 {C 2 \\(l - + C 3 [|(l - W^Wl*} 



We give a similar a priori error estimate for the BQNL method in Theorem 15.21 Following the 
notation adopted in Theorem 15.11 let ^> t ° tal be the energy given by ^ tal (y) := ^ q ^ l (y)— < f,y >. 
Let y be a minimizer of the energy <& total such that A := ming A^ > for 5 2 §(y). 
Theorem 5.2 (A priori error estimate for the BQNL method). Let <3?^°* a ' ; y ; and A be as above. 
Assume that A > 0, and that min^y^ > ^- > 0. There exist constants Si := o~i(A) and 5 2 ■= 
S 2 (mm^y'^, A,C 2 ,C 3 , /3) so that if 

2eC 3 \\A^y'l\\^ + 2 e 2 {C 3 ||(l - ft_i)j/f \\t» + C 4 ||(l - ^)(y'l) 2 \\e^} < S x , and 
elc 2 \\A/3^\\ e 2 + e§{C 2 ||(l - fr-x)i/{'\\% + C 3 \\(l - ftXi^) 2 ^} < S 2 , 
then there exists a locally unique minimizer yp of ^ t ° tal which satisfies 

\\y ~ VpWui,* < j [eC 2 \\A^y'{hi + e 2 {C 2 \\(l - Pt-i)y'{%2 + C 3 \\(l - ^)(^') 2 |l/?}] • 

Proof. The proof of Theorem 15.21 is similar to the the proof of Theorem I5.lt we simply use the 
modeling and stability estimates for BQNL in place of those for BQCE. □ 

6. Conclusion 

We have proposed a smoothly blended version of the quasicontinuum energy (QCE) which we 
call the blended quasicontinuum energy (BQCE). BQCE blends the atomistic and corresponding 
Cauchy-Born continuum energies used in the QCE method over an interfacial region whose thickness 
is a small number k of blended atoms. We analyze the accuracy of BQCE applied to the problem of a 
one-dimensional chain with next-nearest neighbor interactions. For this test problem, we show how 
to choose the optimal blending function for weighting the atomistic and Cauchy-Born continuum 
energies. If BQCE is implemented using the optimal blending function, the critical strain error of 
BQCE can be reduced by a factor of k 2 . Thus, we believe that BQCE could be used to accurately 
compute the deformation of crystalline solids up to lattice instabilities. 

We are continuing the development of the BQCE energy by modifying the code developed to 
study the accuracy of quasicontinuum methods for two benchmark problems — the stability of a 
Lomer dislocation pair under shear and the stability of a lattice to plastic slip under tensile load- 
ing |40| . We note that the potential to significantly improve the accuracy of existing quasicontinuum 
codes by easily implemented modifications is a very desirable feature of the BQCE approach. We 
think our theoretical analysis of the accuracy near instabilities for one-dimensional model problems 
can successfully explain most of the computational results for these multi-dimensional benchmark 
problems. 

We expect that some form of the stability and error estimates in this paper can be generalized 
to atomistic models with finite range interactions by extending the techniques given in |22| . We are 
investigating the extension of our one-dimensional analysis to the multi-dimensional setting, but 
we expect that any multi-dimensional analysis would likely be restricted to perturbations from the 
ground state which are far from lattice instability. We will thus need to rely on our one-dimensional 
analysis to attempt to understand the computational results from our multi-dimensional benchmark 
studies [3D]. 
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